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Abstract 

A model for the simulation of Kirkendall voiding in metallic materials is presented based on vacancy diffusion, 
elastic-plastic and rate-dependent deformation of spherical voids. Starting with a phenomenological explanation of 
the Kirkendall effect we briefly discuss the consequences on the reliability of microelectronics. A constitutive model 
for void nucleation and growth is presented, which can be used to predict the temporal development of voids in solder 
joints during thermal cycling. We end with exemplary numerical studies and discuss the potential of the results for 
the failure analysis of joining connections. 
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"Get your facts first, 
and then you can distort them." 
Mark Twain, (1835-1910) 



1. Introduction 



Microelectronic circuit units consist of the func- 
tional chip unit itself and its packaging, which in- 
cludes several electro- mechanical connections, e.g., 
solder joints between different metal layers. Failure 
of these metallic components is a well established 
cause of failure of the whole microelectronic system. 
The general setup of a typical microelectronic sys- 
tem is illustrated in Figure 1 exemplarily for the 
Flip Chip packaging. Solder balls as well as small. 
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nowadays lead-free joints which are typically made 
of Sn-containing alloys {e.g., Sn-Ag or Sn-Ag-Cu) 
hold the multi-layered unit in position. In addition, 
the solder joints provide electrical conductivity be- 
tween the coppered layers. "Aging" of the solder al- 
loy, such as phase separation, coarsening or the for- 
mation of InterMetallic Compounds (IMCs), as well 
as the formation and growth of pores and cracks in 
the vicinity of heterogeneities significantly effect the 
life expectation of the joints and considerably influ- 
ence the reliability of the whole component. 

During manufacturing the (molten) Sn-rich solder 
wets the copper pad and IMCs are formed due to an 
interfacial reaction [8j. In particular, the copper rich 
CusSn is expected to grow adjacent to the copper 
substrate and CueSns will form adjacent to the Sn- 
bascd solder, c/., Figurc[21 However, in practise both 
IMCs form together with further intcrmctallics to 
irregular shaped layers of initially 2-5 /im hight. Due 
to the rcflow process and to the thermal cycling dur- 
ing service the solder joint ages and the IMCs grow 
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Fig. 1. Intermetallic Compounds at copper-solder interfaces 
in microelectronic components, micrograph courtesy of K. 
Miiller, Neue Materialien Bayreuth, Germany, 2005. 

and may reach a thickness of 20 fim and more, c/., 
[17lll| . Another consequence of solder jomt aging is 
the condensation and growth of so-called Kirkendall 
voids (mainly) within the intermetallic zones. The 
physics behind this mechanism may be sketched as 
follows: neighboring phases or compounds change 
in a way that the volume of one region grows and 
the volume of the another phase reduces. In case of 
Sn-based solders such regions arc typically the in- 
termetallic compounds CusSn and CugSns. These 
IMCs show different diffusion coefficients w.r.t. Cu 
and, therefore, the diffusion of Cu from the pad via 
the interface Cu/CusSn into CusSn is much slower 
than the diffusion of Cu from CuaSn into the CugSns 
scallops, which also cannot be "corrected" by the 
invers diffusion of Sn through the CugSns/Cus in- 
terface. Because of the unbalanced Cu-Sn diffusion 
vacancy-sized voids are left which coalesce to Kirk- 
endall voids, c/.. Figure m Additional vacancies and 
defects in the crystal lattices are generated by plastic 
deformation of the solder material and assist in the 
process of void growth and material degradation. 

The scope of this paper is to model the conden- 
sation and growth of such voids in IMCs within a 
continuum mechanical framework. To this end we 
idealize the material as a homogenous medium with 
arbitrarily distributed vacancies and study the for- 
mation of voids and their growth up to significant 
size. In the following section wc outline the consti- 
tutive equations which enables us to model the void 
condensation and growth process. Subsequently, we 
study short term and long term effects and present 
exemplary results. Finally the temporal develop- 
ment of an ensemble of differently sized voids is 
investigated by means of a void size distribution 
function. 




Fig. 2. First row: Kirkendall voids witliin tlie CusSn layer, 
pliotograpli from |10) . Second row: Crack initiation by void 
coalescence after 1000 thermal cycles between -40° C and 
125° C, photograph from [17]. 

2. Void nucleation by vacancy diffusion 
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Fig. 3. Model of vacancy diffusion and model of a single void 
before and after deformation. 

To explore the feasibility of void condensation 
out of the scale of crystal lattice defects as a void- 
nucleating mechanism wc employ here a model of 
vacancy condensation. To this end let us consider 
a small void possibly just a few atomic spacings in 
diameter and surrounded by a supersaturated back- 
ground vacancy concentration Cb generated, e.g., by 
unbalanced diffusion in the sense of Kirkendall. To 
render the problem analytically tractable, we ide- 
alize the defect to be spherical with characteristic 
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radius a, the diffusion constants to be isotropic, and 
we assume that a steady state vacancy concentra- 
tion profile is maintained all times, i.e., dcyac/dt <C 
dJvac/dx, where Jvac is the vacancy flux. This 
eliminates the time dependence of the solution and 
confers spherical symmetry to the problem. With b 
being the (large) radius of the basin of attraction 
around the void, c/.. Figure [3l and b/a — > oo, the 
diffusion equation for volume concentration c re- 
duces to the following boundary value problem in 
spherical coordinates: 



dr V dr / 
subject to the boundary conditions 

c(r = fe) = Cfc , 

c(r = a) = Ca = Cqc'^I"- . 



(1) 

(2) 
(3) 



Here the equilibrium vacancy concentration near a 
free surface is given by 



Co 



(4) 



and the concentration at the void surface follows 



as above with d ■ 



2-1 Vy 

kT 



due to the Gibbs-Thomson 



effect. In these expressions 7 is the surface tension, 
Vv the atomic volume, k the Bolzmann constant and 
Ey is the free-energy gain/loss resulting from adding 
a vacancy into the system. The solution of ([I])- ([3]) is 
elementary, namely 

c(r) = Cb- (cb - Ca)- . (5) 
r 

For void growth to take place there must be a net 
flux of vacancies into the void, which requires Cb > 
Ca- This in turn requires 



/27W 
Cb > coexpl — — 

V anfci 



(6) 



For very small values of ag or Cb this inequality is 
not satisfied and voids fail to grow. However, the 
value of radius oq which equals relation ([6]) at a given 
value of Cb marks the inception of void growth, i.e., 
a critical nucleation size. Furthermore, the flux J is 
defined to be the change of volume per unit area and 
time, viz. 



J{r) = 



1 



A-KO^ dt 



—Tia" 
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(7) 



dc 



Applying additionally Fick's law J(r) = —Dv-q;^ 
to Eq. ([7]) yields for voids of radius a by means of 
Eq. ®: 

Dv t _ _ _d/a\ 



a{a) 



[Cb - CqC 



(8) 



where Dy is the vacancy diffusion coefficient. As 
the void grows, Ca decreases according to Eq. ([3]) 
which amplifies the concentration gradient towards 
the void and, thus, accelerates the void growth. The 
rate of void growth may be computed by means of 
Eq. jll). 

Following the strategy of Wagner [T3] but allowing 
for an additional source term s{t) the background 
vacancy concentration Cb can be subjected to a "void 
volume balance" of the form: 



d{a, t) d{a) Ana^ da 



= s{t)- / d{a,t) Dv{cb - coe'^/'')ATrada. (9) 



Here d{a, t) denotes a void size distribution function 
describing the fraction of voids with a specific size 
a G [avac, Omax] ^t time t. At this point we neglect the 
statistics and consider only one void size, d{a, t) = 
1, but we will study different void distributions in 
Section[ni The source term s{t) in Eq. ^ represents 
a vacancy production rate due to unbalanced diffu- 
sion which is caused by the different diffusion coeffi- 
cients of the migrating substances. Please note that 
Eq. ([8]) has been derived in [5] in a thermodynamic 
consistent manner (for a constant number of vacan- 
cies). In particular, the vacancy diffusion coefficient 
Dv has been identified as the tracer diffusion coef- 
ficient of vacancies. 

In addition we ask for an energetic formulation 
of the above vacancy diffusion problem. In particu- 
lar, we look for a diffusion rate potential <I'(a, a), for 
which the variational form: 



(5$(d, a) 
5a 







(10) 



holds. For that reason we follow the strategy in [4] 
and multiply Eq. ([8]) with a characteristic factor 
Ey/Dy. A subsequent integration w.r.t. d finally re- 
sults in 



$((2, a) 
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E„d , 



I Cfc - CqC 



d/a\ 



(11) 



3. Constitutive model of void growth 



Once voids are nucleated diffusion is not the only 
mechanism which triggers their growth within the 
IMCs. To set up a general variational model for 
void growth in a deforming material we postulate 
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the existence of a free energy density function as- 
sociated with the deformation of expanding voids 
and embedding material. AdditionaUy we require 
the time-dependent constitutive equations to derive 
from power potentials. Thence we assume the power 
of the external forces acting on the material to be 
completely compensated by the change of its free en- 
ergy and its rate potentials. 

Let us now ask, which energy contributions result 
from the deformation of one void subjected to (a), 
the power P of a remote pressure p(t). These are: 
(b) the energy of the free void surface 5", (c) the de- 
formation energy W of the embedding material and 
(d) the rate power of creep deformation ^! and diffu- 
sion $. Additional energy and power contributions 
may play a role in specific regimes, e.g., the kinetic 
energy in case of a very rapid loading, c/., |16] . Now 
an action integral I{a) can be formulated as sum of 
all rate of energy and power contributions. Hamil- 
ton's principle simply requires stationarity of the 
action integral 51{a) = 0, or, equivalently, 

5_ 
Sd 

In what follows we imagine the material to be a 
conglomerate of (initially very small) spherical voids 
each at every instance completely embedded in the 
material, i.e., we exclude here the process of coales- 
cence, c/., FigureO Consider now a void surrounded 
by a sphere of influence of material with radius b and 
let it expand for some reason. Presuming a volume 
preserving deformation it holds for all r e [a, b] 



dt ^ 



5) + * + $ - P 



= 0. 



(12) 



d 47r^ 3 

(r 

dt 3 ^ 



(13) 



in which the index refers to the initial state, c/.. 
Figure [3l The kinematic relation p3)) will be em- 
ployed subsequently to express functions of radius b 
as functions of current void radius a and the initial 
geometry. Furthermore, the rate of straining of the 
void surrounding material can be defined as 



(14) 



. (dcf) dr d (a .\ 2a . 
dr dr \r^ J r^ ' 
where r is the radius of void surrounding material, 
see Figure El 

In what follows we consider the different energy 
contributions (a-d) in detail. 

3.1. External power 

The external power put into the system by an 
applied (positive or negative) pressure p{t) reads 



P^o. = ^J^Pit)dV 



and for one void we obtain 

P{d,a) = p{t)ATra'^d. (15) 

In the solder material under consideration here the 
effect of thermal cycling on the growth of voids is of 
particular interest. To this end temperature cycles 
T{t) are considered to strain the material and to 
induce a pressure of the form 

pit) = K*{a)3a[To~T{t)] (16) 

where k* is the effective bulk modulus and a is the 
thermal expansion coefficient. Following a classical 
approach of homogenization, cf., [917] , we obtain 
for an assemblage of spherical voids enclosed in an 
isotropic material with bulk modulus k and shear 
modulus fi the effective bulk modulus k* as 

2^ 3k 4- 4/2 



1 



63 Sn + Afib^/a^ 



(17) 



The temperature cycles between — 40°C and 
125°C within one hour, in which 15 minutes at up 
soak and low soak temperature, respectively, see 
Figure HI 
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Fig. 4. Thermal cycling between -40° C and 125° C, cycles of 
60 min with 15 min hold time. 



3.2. Surface energy 

The surface energy of one void with radius a{t) is 
written as 

S{a) ^ Ana^-f , (18) 

where 7 is the surface energy per unit undeformed 
area [N/m]. 

3.3. Elastic-plastic deformation energy 

The deformation energy for the elastic-plastic ma- 
terial response can be derived as follows. We pre- 
sume a Ramberg-Osgood power law, 

where a and e are the effective stress and strain, re- 
spectively. Furthermore e'^ denotes the elastic strain 
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component andE = E{T) and do = <^o{T) represent 
temperature dependent Young's modulus and ini- 
tial yield stress. Exponent n S [1, oo) determines the 
stress-strain curve; in particular n — I prescribes 
linear elasticity, whereas n — > cx) enforces perfect 
(rigid) plasticity. In order to resolve Eq. (fT9|) let us 
assume the elastic strain to be given and decompose 
e = e'^ + €P. Then the dissipated energy of the defor- 
mation per unit volume can be computed from 



ayeP dt 



(7y deP 



nao f E 



n + 1 \ (To 



(20) 



and with reference strain eo = (Jo/ E the dissipated 
deformation energy for one void with surrounding 
material follows as 



W{a-T) 



a n+1 



47rrMr . (21) 



Let us now refer to the kinematic relation in Eq. (|14p 
and assume that the void radius history a(t) grows 
monotonically from ao to ai, then decreases mono- 
tonically from ai to 02, and so on. The integration 
of Eq. p4)l with respect to the time gives 



e(^o,i) = g log 
2, 

+ o log 



■log 





f r3 


-al 




rl 
Vrl 




0.2'^ - 
03^ - 


Vrl 


-al 


- 




-al 



(22) 



and grouping terms corresponding to increasing and 
decreasing intervals yields 



^{ro,t) = - log 



m + '^o' - ao 



3 ^\q^{t) + rl-al, 



(23) 



where q{t) is the maximum radius attained by voids 
of the current size a. 



qit) 



max a(r) 

0<r<t 



(24) 



i.e., q{t) is a monotonically increasing function for 
every history of a(t). Then expression in Eq. (|23|1 
can be summarized, viz. 



e('"o,i) = 3 log 



{q^it)+rl-al) 
rl {a^t) + rl-al) 



2 1 



(25) 



and with Eq. (|13p the dissipated deformation energy 
in a shell surrounding the void is 

W{a,q;T) = in 



1 



n 



2 

3eo 



X log 



(qitf 



rl{a{t)3 + rl~al) 



eo 



Tq dr 



(26) 



It is worth mentioning that the function in Eq. (|26p 
is trackable by analytical means only in the special 
case of n — > 00. Moreover, in moderate hardening 
materials we know the elastic strain component to 
be e*^ « eo. Therefore we set the last term in brackets 
to one but guarantee the energy to be not negative, 
i.e., W >0. 

3.4. Rate- dependent deformation (Creep) 

The creep-power potential is formulated in order 
to capture the effect of the rate of deformation on 
the rate of straining. Experimental observations re- 
ported, e.g., in [1], indicate that the material near a 
void is subjected to a state of stress that is likely to 
cause power-law creep, i.e., 

^) e,exp(-Q,/i?T) , (27) 

where e is the strain rate, m is a creep exponent 
and e'c is a material Constance. The temperature 
dependence of the strain rate is controlled by the 
thermal activation energy Q^. Here we summarize 
the last terms in Eq. (|27p to a reference strain rate, 
eo — e'c cxp(—Qc/RT), with small values of eo corre- 
sponding to creep dominated deformation and eo 
00 to a time independent behavior. The creep-power 
potential per unit volume is defined by 

mcjoe / e 



, d e = 



1 



1 



(28) 



For simplicity we assume now a linear rate depen- 
dence, m = 1, and by integration over the volume 
follows the creep potential for one spherical shell as 

*(d, a) = / -PI- ( J"' 1 47rr^dr 
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(29) 
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where we again make use of Eq. p4p . 

Setting the external power equal the "sum of the 
internal powers" according to ansatz p2|) yields an 
ordinary differential equation for the void size a{t). 
In particular, for growing voids and for n — > oo in 
Eq. ([26|) we obtain the expression 



= - ^■K-p{t)a^ + — (Tofl^ log 



3eo 
87ra7 



1 

7^ (^f- 



■ Ca) 



(30) 



which can be solved (numerically) for all different 
void sizes of interest. We will outline selected results 
subsequently. 

4. Nucleation of voids out of small defects 

At first we study the nucleation of voids, i.e., the 
formation of pores out of vacancy sized defects. It 
can easily be seen from Eq. pO[) that for very small 
values of ao the diffusion term dominates. In the ini- 
tial stages of void growth the elastic-plastic material 
response is of minor influence. As well we neglect 
the external loading for a start and the equation of 
motion reads: 



(Cfc 

a 



(31) 



Thus the void condensation process is driven by an 
interplay of vacancy diffusion and surface energy 
contributions. In particular, for Eq. ([6]) to be satis- 
fied the initial size of the defect must exceed a crit- 
ical nucleation radius which is here approximately 
2-3 times the vacancy size (0.25 nm). To keep voids 
growing a vacancy production source is required, 
otherwise the voids either collapse or reach a steady 
state — depending on the magnitude of surface ten- 
sion. In the material under consideration here the 
vacancy production results from the unbalanced dif- 
fusion, i.e., the different diffusion coefficients of the 
migrating substances, [2]. In a first approximation 
its rate is assumed to be proportional to the differ- 
ence of the diffusion coefficients: 



s{t) ^ CK{t){Dc^ 



D 



Sn) 



(32) 



where -Dcu and Dsn are the diffusion coefficients of 
copper into tin and vice versa, ck is a weighting 
factor with unit [m^] which accounts for progressing 
IMC growth. Note that plastic deformation is an 
additional source of vacancy production and factor 



Ck may also depend on time via the rate of plastic 
straining. 
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Fig. 5. Condensation and growth of voids with different ini- 
tial radii ao driven by an interplay of diffusion and surface 
energy contributions. 

Table 1 

Material constants for diffusion 





CO 




d [m] 


s [1/s] 


7 [N/m] 


10-" 


io-'5 


10-* 


5- 10-'5 


10-7 


1 



Although Cu-Sn intermetallics are of great inter- 
est for the electronic packaging industry the avail- 
able material data on such compounds vary consid- 
erably, see, e.g., |3I6I13| and references therein. With 
typical values summarized in Table [T] we obtain for 
the solution of Eq. ([31]) the results'^ displayed in 
Figure[5l Small defects with an initial radius of ao = 
0.5 nm (twice the size of a vacancy in copper) col- 
lapse immediately. Defects of size greater than nu- 
cleation size will grow with a rate of void growth 
depending on the vacancy production rate. On the 
other hand, big voids, e.^., ao = 500 nm, basically 
fail to grow, here the influence of diffusion is to small 
for significant void growth. 

Figure [6] shows the same void growth model but 
now with external loading, i.e., the voids are sub- 
jected to a pressure history induced by thermal cy- 
cling. On relatively small voids this pressure has lit- 
tle effect, their growth is in first instance given by 
diffusional effects. However, if the void size reaches a 
significant magnitude the void grows unbounded - a 
clearly non-physical effect - which shows the neces- 
sity to incorporate elastic-plastic deformation into 
the model. This effect should be considered in the 
following section. 



The ODEs of Eq. (I31I I have been solved numerically using 
the Matlab® solver 0DE45 as well as an explicit 2th order 
RUNGE-KUTTA procedure. 
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Fig. 6. Influence of the external load induced by subjecting 
the material to the temperature cycles of Figure |4] on the 
diffusion driven void growth of two voids. 

5. Growth of voids due to thermal cycling 

The full model of Section [5] is now studied for 
medium sized and big voids with material param- 
eter given in Tabled Note that the elastic-plastic 
deformation energy contribution in Eq. (j30[) is sim- 
plified in order to get an analytical expression; a 
computation of the full model requires a numerical 
integration of Eq. ([26|) . In Figure[7|we see the differ- 
ent rates of growth for medium sized and big voids. 
Here we neglect creep effects at first, i.e., eo — > oo. 

Both void sizes are subjected to the same thermal 

Table 2 



Elastic-plastic material constants * 



E [GPa] 


K [GPa] 


At [GPa] 


(JO [MPa] 


a [1/K] 


100 


80 


50 


450 


19 • lO"** 



cycling of Figure 01 but their evolution history de- 
pends strongly on the initial void size. In particular, 
smaller voids clearly grow slower than their initially 
big companions. That supports the experimental ob- 
servation of several (relatively) large Kirkendall 
voids in IlVICs. The differences in the void size vs. 
time response for the two cases of elastic-plastic ma- 
terial behavior (n = 5) and perfect plastic approxi- 
mation [n — > oo) are small, the final void size is ba- 
sically determined by the external loading. There- 
fore, a simplified approximation as given in Eq. pO|) 
seems to be justified. Let us point out that the dif- 
fusion effects invoke a "smoothening" of the void 
growth curves. Even if the underlying elastic-plastic 
theory is time independent, the material response 
has no sharp edges as the load history would suggest. 

* Experimental data courtesy of Prof. Miiller, Inst. f. Mech. 
(LKM), TU Berlin, Germany 
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Fig. 7. Growth of a medium sized and large void induced by 
thermal cycling for two elastic-plastic material laws accord- 
ing to Eq. 1191 n = 5 and perfect plasticity. 

Figure [8] displays the void size vs. time for differ- 
ent reference strain rates. Small values of eo clearly 
damp the evolution of voids. However, after enough 
time has passed the final void size will reach the 
same value as in the time-independent case as can 
be seen from the curves for eo = 10"^ and eo = 10"^ 
in Figure [51 
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Fig. 8. Effect of creep on the growth of voids of size 
ao = 500nm for different reference strain rates. 

6. Evolution of a Void Ensemble 

In order to investigate an ensemble of voids with 
different initial radii we introduce ~ as already in- 
dicated in Section [5]- a void size distribution func- 
tion d{a,x,t). To derive an evolution equation for 
the void distribution we make use of a mesoscopic 
concept described in detail, e.g., in [12115116] . In 
this work we do not account for a dependence of 
the void distribution on the spacial position x, and 
set d(a,x,t) = d{a,t) = d. Then we establish for a 
constant number of voids a balance equation for the 
void size distribution d in the following form: 
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dd d ~ d , . , s 







(33) 



The value of a(a) can be calculated from the consti- 
tutive model introduced in Section [2] and [3l More- 
over we neglect the spatial velocity of the material 
element under consideration, w 0, so that the sec- 
ond term in Eq. (j33p vanishes. 
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Fig. 9. Comparison of the numerical solution of Eq. I I33I I and 
I I34I I for 5 different radii after 2 thermal cycles. 

To solve Eq. ([55]) numerically, the spatial and tem- 
poral derivatives were discretized by finite Differ- 
ences. The solution of Eq. (|30[) is obtained by means 
of a 2nd order Runge-Kutta method. For ease of 
computation we also solved the quasi-linear pendant 
of Eq. dSSl), viz. 



dd . , .dd 



(34) 



As is shown exemplarily in Figure [9] both equations 
yield almost identical results. However, to solve the 
full Eq. (|33p a significant finer spatial discretization 
is required. Therefore we proceed here with solutions 
of Eq. Furthermore, for a better understanding 
of the results we neglect any void production, i.e., 
s{t) =0 in Eq. ([9]). However, there is no extra effort 
to include the production term into the simulations 
and in Eqs. and (pi]) . 

We start with a discretized GAUSSian distribu- 
tion of 15 different void radii, a = 10 . . . 150 nm. 
Within several stages of thermal cycling we obtain 
the temporal development of d as illustrated in Fig- 
ure llOl The initially symmetric distribution changes 
to an asymmetrical distribution in such a way that 
the fraction of smaller voids decreases and the big- 
ger voids grow. Such results are well known from 
so-called LSW-theories for Ostwald ripening, [14], 
where bigger "grains" grow at the expense of the 
smaller ones due to the Gibbs-Thomson effect. In 
our model this effect also appears in the first stages 
of void growth, in which vacancy diffusion domi- 
nates. During proceeding growth the voids reach a 
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Fig. 10. Development of an initial GAUSSian void distribution 
after 1, 2, 5 and 10 thermal cycles. 



size, for which the evolution is characterized primar- 
ily by elastic-plastic deformation. For such stages 
the distribution function considerably differs from 
a typical LSW distribution. In particular the num- 
ber of large voids extremely exceeds the number of 
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small voids. 
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Fig. 11. 1st row: short time behavior of d. 2nd row: long 
time behavior of d. 

A different view on the evolution of void size dis- 
tribution shows Figure [TT] The first plot displays 
the size distribution of small voids in the beginning 
of thermal loading. In the long range (second plot in 
Figure [TT|) the growth of the bigger voids driven by 
clastic-plastic deformation by far exceeds the growth 
of small voids. The final distribution substantially 
differs from the initial GAUSSian form. In particular 
the void regimes primarily consists of large voids, 
which may allow the assumption of crack formation 
by void coalescence in the immediate future (c/., 
Figure [21 2nd row). 

At second we study a distribution of void sizes, 
d, which - in a more realistic manner - accounts for 
the dominance of small voids in the initial state. To 
this end we assume an exponential distribution of 25 
voids with a = 10 . . . 250 nm. As displayed in Fig- 
ure [Inland [12] we observe again a growth of the frac- 
tion of big voids on cost of the smaller ones. In par- 
ticular, the void size a reaches large values and the 
distribution function is stretched over a wide range 
of void sizes. Here wc stopped the simulations after 
a void volume fraction of « | is reached. Note that 
our model does not account for void growth due to 
coalescence. Unfortunately this limits the predictive 
capabilities of the void size distribution analysis. 



However, already a significant amount of large voids 
in a material clearly indicates the onset of failure. 
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Fig. 12. Development of an initially exponentially decreasing 
void size distribution after 1, 2, 5 and 10 thermal cycles. 
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Fig. 13. Compact illustration of the evolution of the distri- 
bution function for different numbers of thermal cycles. 



7. Conclusions 

We presented a constitutive model to predict the 
condensation and growth of Kirkendall voids in 
elastic-plastic metals, in particular in the IMCs oc- 
curring at the interface of microelectronic solder 
joints. To this end the influence of vacancy diffusion, 
surface energy, elastic-plastic and creep deformation 
on the evolution of void ensembles during thermal 
cycling was investigated. It turns out that nano- 
voids collapse, whereas voids which are small but ex- 
ceed a critical radius (of a few vacancies) grow driven 
by diffusional effects. On the other hand the growth 
of bigger voids is primarily driven by elastic-plastic 
deformation of the void surrounding material. We 
found that work hardening plays a minor role and, 
as expected, creep decelerates the void growth. 

Furthermore we studied the temporal develop- 
ment of void ensembles under thermal cycling. Here 
the presented model for void growth is employed to 
(numerically) solve the balance of the void size dis- 
tribution function. We found that the evolving dis- 
tribution initially resembles the ones known from 
LSW theories, whereas the distribution during pro- 
ceeded loading evolves such that the amount of large 
voids drastically increases. Such behavior correlates 
to experimental studies on Kirkendall voids. 

The presented results are applicable to derive the 
temporal evolution of the effective properties of in- 
termetallics. The constitutive model shall be incor- 
porated in Finite-Element Analysis tools on the ma- 
terial point level to study the mechanical behavior 
of IMCs under realistic loading regimes. Such ana- 
lyzes then allow for a better prediction of live time 
and strength of microelectronic joining connections. 
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